Similarity submodules and root systems in four dimensions 



Michael BAAKE a [] and Robert V. MooDY b 

a) Institut fiir Theoretische Physik, Universitat Tubingen, 
Auf der Morgenstelle 14, D-72076 Tubingen, Germany 
b) Department of Mathematical Sciences, University of Alberta, 
Edmonton, Alberta T6G 2G1, Canada 



Dedicated to H. S. M. Coxeter 



Abstract 

Lattices and Z-modules in Euclidean space possess an infinitude of subsets 
that are images of the original set under similarity transformation. We classify 
such self-similar images according to their indices for certain 4D examples that 
are related to 4D root systems, both crystallographic and non-crystallographic. 
We encapsulate their statistics in terms of Dirichlet series generating functions 
and derive some of their asymptotic properties. 

Introduction 

This paper begins with the problem of determining the self-similar images of certain 
lattices and Z-modules in four dimensions and ends in the enchanting garden of Cox- 
eter groups, the arithmetic of several quaternion rings, and the asymptotics of their 
associated zeta functions. The main results appear in Theorems || and [3|. 

The symmetries of crystals are of fundamental physical importance and, along 
with the symmetries of lattices, have been studied by mathematicians, crystallogra- 
phers and physicists for ages. The recent interest in quasicrystals, which are non- 
crystallographic yet still highly ordered structures, has naturally led to speculation 
about the role of symmetry in this new context. Here, however, it is apparent that 
a different set of symmetry concepts is appropriate, notably because translational 
symmetry is either entirely lacking or at least considerably restricted in scope. 
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One of the most obvious features of quasicrystals is their tendency to have copious 
inflationary self-similarity. Thus instead of groups of isometries, one is led to semi- 
groups of self-similarities that map a given (infinite) point set into an inflated copy 
lying within itself. Ordinary point symmetries then show up as a small part of this, 
namely as the "units" , i.e. as the (maximal subgroups of) invertible elements. 

The importance of self-similarities is well-known, and has also been used to gain 
insight into the colour symmetries of crystals [[3^, |33| and, more recently, of quasicrys- 
tals 0, |5|, |25[ . The focus in the latter cases was on highly symmetric examples in the 
plane and in 3-space since they are obviously of greatest physical importance. 

In this contribution, we extend the investigation of self-similarities to certain ex- 
ceptional examples in 4-space, namely the hypercubic lattices (of which there are two, 
represented by the primitive hypercubic lattice Z 4 and the root lattice D 4 or, equiv- 
alently, its weight lattice Dl) and the icosian ring, seen as the Z-span of the root 
system of the non-crystallographic Coxeter group if 4 , see JlO], [19], [L3| for notation and 
background material. The case of lattices has recently also been investigated by Con- 
way, Rains and Sloane M who are generally interested in the question under which 
conditions similarity sublattices of a given index exist. Their methods are comple- 
mentary to ours, and more general, but do not seem to give direct access to the full 
combinatorial problem which we can solve here. 

It is useful to digress briefly to discuss the role of H4 in the context of aperiodic 
order. It is a remarkable fact that the non-crystallographic point symmetries relevant 
to essentially all known physical quasicrystals are actually Coxeter groups, namely 
the dihedral groups h^k), k = 5,8,10,12 and the icosahedral group if 3 Apart 
from the remaining infinite series of dihedral groups, the only indecomposable non- 
crystallographic Coxeter group is H± of order 14400, which is the most interesting 
of them all. In spite of being four-dimensional in nature, there are several good 
physical reasons to explore the symmetry expressed by this group, and because of its 
connections with exceptional objects in mathematics, including the root lattice Eg, 
there are good mathematical reasons, too. 

From the point of view of quasicrystals, may be viewed as the top member of 
the series H 2 : = -^(5) C H 3 C of which the first two have been the subject of great 
attention, see |2B| and references therein, while H4 appears as symmetry group of the 
Elser-Sloane quasicrystal [|l7j]. Mathematically, this family belongs together and H4 
is the natural parent of the others. 

Now, if 4 has a natural quaternionic interpretation which arises as follows. The 
group of norm 1 units of the real quaternion algebra is easily identifiable with SU(2), 
see |22] for background material and notation. Using the 2- fold cover of SO (3) by 
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SU(2) we can find (in many ways) a 2-fold cover of the icosahedral group inside SU(2). 
This is the binary icosahedral group / of order 120 [14, p. 69]. The point set I is a 
beautiful object, namely the set of vertices of the exceptional regular 4-polytope called 
the 600-celig |]| Ch. 22] (under the standard topology of IR 4 carried by the quaternion 
algebra). The Z-span of I is a ring I, dubbed by Conway the ring of icosians. This 
ring is closed under complex conjugation and under left and right multiplications by 
elements of /. The group of symmetries of I so obtained is isomorphic to H4 acting as 
a reflection group in IR 4 . The ring I is itself quite a remarkable object. It is naturally 
a rank 4-module over Z[r], r := (1 + \/5)/2, and a rank 8-module over Z (so it is 
certainly dense in the ambient space IR 4 ). In fact, as an aside, it has a canonical 
interpretation as the root lattice of type E 8 with J U rl making up the 240 roots of 
E 8 . Restricting to the pure quaternions puts us in the 3-dimensional icosahedral case, 
and by further restriction we can get the H2 situation. 

Now we can state our problem for the icosian case. We have pointed out that the 
additve group I has a large finite group of rotational symmetries coming from the left 
and right multiplications by elements of I - namely 120 2 /2 = 7200 such symmetries. 
But we have seen that in the study of quasicrystals we have to pay attention to 
self-similarities, too. So we are now also interested in rotation-inflations of IR 4 that 
map I into, but not necessarily onto, itself. Each such self-similarity maps I onto 
some submodule of finite index, and our question is to determine these images and to 
count the number of different similarity submodules of a given index. This leads us to 
introduce a suitable Dirichlet series generating function which encodes the counting 
information and its asymptotic properties all at once, and indeed determining its exact 
form is a number theoretical problem that depends heavily on the fact that I can be 
interpreted as a maximal order in the split quaternionic algebra over the quadratic 
field Q(VE). 

The other situation that we wish to discuss in this paper is crystallographic in 
origin, but it nevertheless is very much the same problem. It is well known that 
there are two hypercubic lattices in 4-space [|3^, |7|, namely the primitive and the 
centred one (face-centred and body-centred are equivalent in 4-space by a similarity 
tranformation). Let us take Z 4 and the root lattice -D4 as suitable representatives. 
Note that they have different holohedries, namely one of order 1152 (denoted by 33/16 
in 0, Fig. 7]) for _D 4 , which coincides with the automorphism group of the underlying 

2 The 600-cell and its dual, the 120-cell, are two of the three exceptional regular polytopes in 
4-space jl^, p. 292]. The remaining one, the 24-cell, also occurs, later in this paper. Beyond 4 
dimensions, the only regular polytopes are the simplices, the hypercubes, and their duals, the cross- 
polytopes (or hyperoctahedra). 
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root system, and an index 3 subgroup (denoted by 32/21 in J7], Fig. 7]) for Z 4 . Due to 
the previous remark, we may take the weight lattice D\ instead of if we wish, and 
we will frequently do so. Given any of these cases, we want to determine the Dirichlet 
series generating function for the sublattices that are self-similar images of it. What 
makes this situation tractable is that, parallel to the icosian case, there is a highly 
structured algebraic and arithmetic object in the background, namely Hurwitz' ring 



of integral quaternions |2(], |TTJ. In our setting, it is J = D\, and it is again a maximal 
order, this time of the quaternionic algebra over the rationals, Q. The results for the 
Hurwitzian and icosian cases are striking in their similarity. Another example is that 
of the maximal order in H(Q(\/2)) the treatment of which equals that of I whence we 
only state the results. 

The structure of the paper is as follows. In the next Section, we set the scene by 
collecting some methods and results from algebra, and algebraic number theory in 
particular. This will be done in slightly greater detail than necessary for a mathemat- 
ical audience, but since there is also considerable interest in this type of problem from 
the physics community, we wish to make the article more self-contained and readable 
this way. The two following Sections give the results, first for lattices, and then for 
modules. We close with a brief discussion of related aspects and provide an Appendix 
with material on the asymptotics of arithmetic functions defined through Dirichlet 
series. 



Preliminaries and Recollections 

We shall need a number of results from algebraic number theory, both commutative 
and non-commutative. First of all, we need, of course, the arithmetic of Z, the ring 
of integers in the field Q. All ideals of Z are principal, and they are of the form 
a = mZ with m G Z. If a ^ 0, the index is [Z : a] = \m\. The corresponding zeta 
function, which can be seen as the Dirichlet series generating function for the number 
of non-zero ideals of a given index, is Riemann's zeta function itself 

1 °° 1 1 

cm = = E-s = Ut^-s (i) 

aCZ 1 J m=l p eV r 

Here, V denotes the set of (rational) primes, and the second representation of the 
Riemann zeta function is its Euler product expansion. It is possible because the 
number of ideals of index m is a multiplicative arithmetic function, a situation that 
we shall encounter throughout the article. 
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Next, we need the analogous objects for the real quadratic field Q(\/5) = Q(t). 
The ring of integers turns out to be 



Z[t] = {m + nr \ m,n G Z} (2) 

where r = (1 + \/5)/2 is the fundamental unit of Z[r], i.e. all units are obtained 
as ±r m with m G Z. Again, Z[r] is a principal ideal domain, and hence a unique 



factorization domain [I8[ ch. 15.4]. The zeta function is the Dedekind zeta function 



36| , §11] defined by 



oo 



acZ[r] 1 1 1 1 m=l 

where o runs through the non-zero ideals of Z[r] and [Z[r] : a] is the norm of a. If 
a = aZ[r], it is given by 

[Z[r] : a] = |N(a)| = \aa'\ (4) 

where ' denotes algebraic conjugation in Q(r), defined by r i— > 1 — r. 
Explicitly, the zeta function reads (see the Appendix for details): 

c Q( r)( s ) = ! _ 5 - s ■ n n _„-^2 ■ n i _ „-2 fl ( 5 ) 

ps±l (5) V F 7 P=±2 (5) F 

= 1 + — + — + — -I 1 1 1 1 1 1 1 1 1 

J- -r 4 « t 5 s -r gs -r lls -r 16 „ -r igs -r 2QS t 25s t 29 s T 31 3 ~ 36 s ~ 41 s ~ 

As before, a(m) is a multiplicative arithmetic function, i.e. a(mn) = a(m)a(n) for 
coprime m, n. It is thus completely specified by its value for m being a prime power, 
and from the Euler product in (|5|) one quickly derives that a(5 r ) = 1 (for r > 0). 
Then, for primes p = ±2 (5), one obtains a(p 2r+1 ) = and a(p 2r ) = 1, while for 
primes p = ±1 (5), the result is a(p r ) = r + 1. 

One benefit of relating the numbers a(m) to zeta functions with well-defined an- 
alytic behaviour is that one can rather easily determine the asymptotic behaviour 
of a(m) from the poles of the zeta function, see the Appendix for a summary. In 
this case, the function a(m) is constant on average, the constant being the residue of 
Cq( t )( s ) a t its right-most pole, s = 1. Explicitly, we get 



1 N 21 ( ) 

j™,]vZX m ) = res s=1 C Q(T) (s) = — |p ^ 0.430409. (6) 

m=l * 
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We shall also need the zeta function of the quadratic field Q(V2), where Z[V2] is 
the corresponding ring of integers, and 1 + y/2 its fundamantal unit [18|]. The zeta 
function reads 

111 

cq(V2)0) = ! _ ■ n n _o-^2 ■ n i _ „-2 fl ( ? ) 

p=±l (8) V F 7 p=±3 (8) F 

= 1 + — + — + — + — + — -I 1 1 1 1 1 1 1 1 

2 s 4 s 7 s 8 s 9 s 14 s 16 s ' 17 s ' 18 s ' 23 s ~ 25 s ' 28 s ' 

The coefficients are a(2 r ) = 1, a(p r ) = r + 1 for p = ±1 (8) and a(p r ) = resp. 
1 for p = ±3 (8) and r odd resp. even. The asymptotic behaviour is given by 
lm%->oo F E m <N <rn) = log(l + y/5)/y/2 ~ 0.623225. 

Let us now move to the non-commutative results we shall need. We will be con- 
cerned with the quaternionic algebra M(K), mainly over the field K = Q or over 
K = Q(t). The case of K = Q(\/2) is treated more as an aside. In all cases, we are 
interested in the corresponding ring of integers, these being the Hurwitzian ring J, the 
icosian ring I, and the cubian ring K. These are maximal orders in their respective 
quaternionic algebras |50, |3~T|. 

Let us first consider J = D\. In terms of the standard basisf] l,i,j,k of H(Q), 
J consists of the points (xo, xi, X2, £3) whose coordinates Xi either all lie in Z or all 
111 Z + \. Though non-commutative, J is still a principal ideal domain, i.e. all left- 
ideals (and also all right-ideals) are principal p0]| . Consequently, we have unique 
factorization up to units, the units being the 24 elements obtained from (±1,0,0,0) 
plus permutations and from |(±1, ±1, ±1, ±1). These 24 units form a group, J x , that 
is isomorphic to the binary tetrahedral group |14], p. 69]. The number of non-zero left- 
ideals (or, equivalently, of non-zero right-ideals) can now be counted explicitly, which 
results in the corresponding zeta function Cj( s ) = SoTFd 7 ' where o runs over the 



non-zero left ideals of J, see |15|, Sec. VII, § 8 and § 9] for details on zeta functions of 



quaternionic algebras. The result is [28 , § 63, A. 15]: 

• The zeta function of Hurwitz' ring of integer quaternions, J, is given by0 

C» = (l-2 1 - 2s )-C(2s)C(2 S -l). (8) 
Using (|l]), one can easily determine the first few terms 

t ( <A — 14_J_4_A4_J_4__6_4_^4__8_4_J_4_13_4_ i JJL _|_ _J_ _|_ (Q\ 

Sjl" 5 ; L 1 43 T gs T lgs T 25s ' 3Qs T 49s T g 4 s 1 S1 3 T 10 qs T 12 is T 144 s T 



The defining relations |11 are: i 2 = j 2 = k 2 = ijk 



4 Note that the formula given in the first line after M Eq. 3.17] contains a misprint in the prefactor. 
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The possible indices of ideals are the squares of integers. It is thus convenient to write 
the Dirichlet series as 

C» = E -hr ( 10 ) 



m 1 

m=l 



so that aj(m) is actually the number of left-ideals of index m (rather than m). This 
then results in cij(2 r ) = 1 (for r > 0) and in a 3 (p r ) = (p r+1 — l)/(p — 1) for odd primes. 
Let us add that a s (m) is also the sum of the odd divisors of m, see |30], sequence M 



3197] or |g sequence A 000593]. 

Let us again briefly comment on the asymptotic behaviour. In this case, the average 
of Oj(m) grows linearly with m, i.e. 

i N 2 

-JjH a ^ m )~Y A N (N^oo) (11) 

m=l 

where the coefficient is half the residue of Cj( s ) a t its right-most pole, s — 1. This can 
easily be calculated from the details provided in the Appendix. 

Note that the linear growth of the average of a,$(m) stems from the definition used 



in (10). With the usual definition (i.e. with denominators m s rather than m 2s ), the 
average would tend to a constant, as in (|6]). 

Now, let us consider the analogous situation, with I being a maximal order in 
the algebra H(Q(r)). Again, all left-ideals (and all right-ideals) of I are principal, 
and we also get unique factorization up to units again |31j . The unit group I x of 
I consists of the 120 elements of the binary icosahedral group / inside I. Taking 
the unit quaternions (1,0,0,0), 1,1,1), |(r, 1, — 1/r, 0) together with all even 
permutations and arbitrary sign flips results in an explicit choice of the group J, and 
hence of I. Again defining the zeta function for one-sided ideals of I we have 

• The zeta function of the icosian ring, I, is 

GOO = ^ (T) (2 a )C Q(r) (2«-l). (12) 
This result follows from |RTL Ch. Ill, Prop. 2.1]. The first few terms of this series read 



- 1 _i_ _i_ _i_ _i_ _2_L_ _i_ _i_ _[1L_ _i_ __ilLL_ _)_ _lLL_ _i_ _ii!_L_ _i_ _li_L_ _i (1%) 

Ml* 5 ; ^ ^ 16 s ~ 25 s ' 81 s ' 121 s T 256 s ^ 361 s ~ 400 s ^ 625 s ^ 841 s ^ 961 s ^ V XJ / 

The possible indices (= denominators) are the squares of integers that are repre- 
sentable by the quadratic form x 2 + xy — y 2 , i.e. of integers all of whose prime factors 
congruent to 2 or 3 (mod 5) occur with even exponent only. Using a definition analo- 
gous to (|10D above, the coefficient a t (m) is again a multiplicative arithmetic function. 
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It is given by a l (5 r ) = (5 r+1 — l)/4 (for r > 0), and, for primes p = ±2 (5), by 
°i(p 2r+1 ) = and a t (p 2r ) = (p 2r+2 — l)/{p 2 — 1). Finally, for p = ±1 (5), one finds 

%(p r ) = E[=o( Z + x )( r - Z + l )P l - lt is now listed as IE sequence A 035282]. 

The asymptotic behaviour of %(m) is similar to that of a,j(m) above, and we obtain 

1 ^ 2 4 1 ( ) 

- %M ~ ^ 375 ^ - 0.249997 -N (N -> oo) (14) 

m=l 

where the slope is again half the residue of £j(s) at its right-most pole, s = 1, see the 
Appendix for details. 

Very similar is the situation of the ring K in H(Q( v / 2)), generated as the Z[V2\- 
span of the basis {1, (1 + i)/y/2, (1 + j)/y/2, + k)/2}. The unit group K x 

is the binary octahedral group of order 48 p. 69], and the symmetry group of K 



contains that of J as an index 2 subgroup. K is again a maximal order and a principal 



ideal domain [|3l|. The zeta function of IK reads [31, Ch. Ill, Prop. 2.1] 



CM =C Q( vs)(2s)Cq(v5)(2s-1) (15) 

_ 1 I 3 I 7 , 16 | 15 i 10 | 48 , 31 , 36 , 30 , 48 , 26 , 
4 s 16 s 49 s 64 s "l" 81 s 196 s 256 s 289 s 324 s ' 529 s 625 s ~* 

and further details can be worked out in complete analogy to the icosian case. 

Let us now briefly describe how the quaternions enter our (mainly geometric) 
picture, and how they provide a parametrization of (S)0(4) = (S)0(4, R), see |22| for 
details. The key is that pairs of quaternions in H(R), i.e. quaternions q = (q , q ± , q 2 , q 3 ) 
as written in the standard basis 1, k of the quaternion algebra over R, induce an 
action on vectors of R 4 via 

M{q 1 ,q 2 )x t = q x xq 2 (16) 

where M(q ± , q 2 ) G Mat(4, R) and x l is x written as a column vector in R 4 . Evidently, 
for nonzero quaternions q 1 , q 2 ,r 1 ,r 2 we have M(q 1 ,q 2 ) = M(r 1 ,r 2 ) if and only if 
r 1 = aq x and r 2 = a~ x q 2 for some a G R\{0}. With q x = (a,b,c,d) and q 2 = 
(t, u, v, w), the matrix M = M(q 1 , q 2 ) reads explicitly 

/ at + bu + cv + dw —bt + au + dv — cw —ct — du + av + bw —dt+cu — bv + aw \ 

bt — au + dv — cw at + bu — cv — dw —dt+cu + bv — aw ct + du+av + bw 

ct — du — av + bw dt+cu + bv + aw at — bu + cv — dw -bt — au + dv+cw 

y dt+cu — bv — aw —ct+du — av + bw bt + au + dv + cw at — bu — cv + dw J 

It has determinant 

det(M( 9l ,qr 2 )) = (a 2 + b 2 + c 2 + d 2 ) 2 ■ (t 2 + u 2 + v 2 + w 2 ) 2 

= {(9i) 2 " (q 2 f} 2 (18) 



(17) 
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and also fulfils 

MM 1 = Vdet M • 1 4 . (19) 

Consequently, when |g 1 | = \q 2 \ = 1, we obtain a 4D rotation matrix and the 
homomorphism M : S 3 x S 3 — > SO (4) provides the standard double cover of the 
rotation group SO(4) [p2]] , with M(q 1 , q 2 ) = M(—q 1 , —q 2 ). The orientation reversing 
transformations, i.e. the elements of 0(4)\SO(4), are obtained by the mapping x \— > 
q 1 xq 2 with unit quaternions q t , q 2 . Let us finally note that, for non-zero quaternions, 

R( qi ,q 2 ) := MA,gr) = -±-M(q v q 2 ) (20) 
|yil \ ( l2\ |yiy2l 

always gives a rotation matrix, which is handy for finding suitable parametrizations of 
groups such as SO (4, Q) or SO (4, Q(r)) in closely related problems, see for details. 



Arguments in common 



In this Section, we focus on I versus J and carry the arguments as far as possible 
without having to separate the two rings I and J too seriously^. Thus we introduce 
the following notation to cover both situations simultaneously: 



K 
o 
O 
C 



z 
J 
z 4 



or Q( v / 5) 
or Z[r] 
or I 

or ZH 4 



(21) 



By L we really mean the ring ol + oi + oj + ok C O, and we observe that 

20 C £ C O. 



(22) 



Let us first note that O is a maximal order in M.(K ) . As such, each prime ideal 
of O corresponds to one prime ideal of o, namely to p := ^3 fl e>, and this sets up a 



one-to-one correspondence between their prime ideals [27], Thm. 22.4]. Furthermore, 
we have unique factorization of each 2-sided ideal 21 of O: 



21 



W ■ . . . • % 



(23) 



where 



fci 



, are prime ideals and all ki > 0. 



5 The case O = K is entirely parallel to that of O 
mention details later when we need them. 



and need not be spelled out here. We shall 
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The prime 2 (which is prime both in Z and in Z[r]) plays a special role in this 
paper. In the case of I, 21 is the prime ideal of I lying over 2Z[r]. The case of J, 
however, is more complicated. Here, (1 + i)J = J(l + i) is the prime ideal lying over 
2Z and (1 + *) 2 J = 2JJ ||20|| . It is this ramification of 2Z in J that accounts for the stray 
factor in Eq. (H) and, later on, in Eq. (|4|). To cope with this, we shall call a G J an 
odd element of J if a ^ (1 + i)J. It is useful to note that a G J is odd if and only if 
\a\ 2 G Z (its quaternionic norm) is odd. 

Let a G O. We define 

A(a) := {r G K | ra G 0} . (24) 

For a ^ 0, o C v4(<x) C (9a, -1 fl i^. Thus -A(a) is a finitely generated o-module and 
hence is a fractional ideal of o |2~I| , Ch. 1.4]. It follows that A(a) _1 is an ordinary ideal 
of o and, consequently, A(a)^ 1 = oc for some c G o. We call c = c(a) G e> (which is 
determined up to a unit of o) the content of a: A(a) = oc(a) . 

Definition 1 We say that a G O is O-primitive if A(a) = o. 

We have just seen that a is (9-primitive if and only if c(a) is a unit in o. So, we have 

Lemma 1 For any a G 0\{O} ; c(a)~ 1 a is a primitive element of O. □ 

Let us now come to the link between submodules and similarities. 

Definition 2 Let L be a 'L-module in M 4 that spans M 4 . M. C L is called a similarity 
submodule (SSM) of L if there is an a & K\{0} and an R G 0(4) so that M. = 
aR(0) C L. If L is a lattice, we call M. a similarity sublattice (SSL) of L. 

Consider the case when M. is an SSM of O. It is immediate that such an M. is an 
o-submodule of O. Since o is a principal ideal domain and O is a free o-module of 
rank 4, we see that M. is also a free o-module of rank 4. Consequently, the index 
[O : M] of M in O is finite. 

We now come to the first crucial assertion in our classification of the similarity 
submodules according to their indices. 

Proposition 1 Let Ai C O be a similarity submodule. Then there exist a, 6 G O, 
with a primitive, such that M. = aOb. In addition, in the case O = J, we can arrange 
for a to be odd. 

Proof: By assumption, M = aR(0) with a e K, a ^ 0, and R G 0(4) = 0(4, R). 
Using ([TED, and noting that O = O, we can write A4 = aOb where a = (a , a 1( a 2 , a 3 ), 
6 = (6 , 6 ls 6 2 , 6 3 ) and a, 6 G M(R). 
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Since 2a{l, i,j, k}b C 2aOb C £, we can infer from the explicit matrix form (|T7|) 
that all the matrix entries of 2M(a, b) lie in o. Combining suitable entries, four at a 
time, we can obtain that 

Sa^eo, for all i, j G {0,1,2,3}. (25) 

Since aOb = arOr~ 1 b for all r G M\{0}, we can arrange that some a { G K\{0}, 



whereupon we see, via (25), that we can choose a,b G M(K) without loss of generality. 
Clearing denominators and using Lemma [I], we may further assume that M. = aOb 
with a G O and a primitive. From (^5]), we now get 8b j a G 0, whence 86^ G A(a) = o 
for all j G {0, 1, 2, 3}. We conclude that 8b G O. 

We now have to dispose of the factor 8 in order to prove the Proposition. Consider 
the 2-sided ideal OaObO of O. Since O is a maximal order 0, f27J [31] in H(K), we 
have a unique factorization 

OaO&O = qj* 1 • . . . • ^ , jfej, . . . , £v G N , (26) 

where , . . . , ^7 are p r i m e ideals of O and D o = p i are in one-to-one corre- 
spondence with distinct prime ideals of o [^7| . Similarly, we have OaO = Y[ anc ^ 
ObO = [] V1\ with ki = mi + m for all i G {1, . . . , r}, m, G N , G Z. 

Since 80bO C (9, the only primes for which rij < is possible are those lying over 
2 G o. Recall that there is exactly one prime ideal of O that corresponds to 2, and 
this is = xO = Ox, where x = (1 + i) in the case (9 = J and x = 2 in the case 
= 1 

Thus we have rii > for all i > 1. If we can now show that n\ — 0, then b £ O 
and our assertion follows. 

Suppose ni < 0. Then mi > |ni| > 0, and we can write 

aOb = ax-^Ox^b (27) 

and Oax^O C Sp^ 1-1 " 1 '^ ..... ^ c C. Similarly, Ox^bO C 0. 

In the icosian case (where x = 2 G e>), the primitivity of a rules out that aaH ni ' G 
0, so ni < is impossible here. In the Hurwitzian case, since (1 + i) 2 = 2 (up to 
units), ni = —1 is still possible. Then, ax l ,xb G and ax -1 is still primitive. We 
take these as the new a, b. This achieves the correct form, M. = aOb, with a,b £ O 
and a (9-primitive. When O = J and a is even, we have a G (1 + i)C?\2C?, and we 
may replace a(96 by ax l Oxb. □ 

Remark 1: Given M. = a 1 Ob 1 C where a 1 ,b 1 G H(iC), the above argument 
shows that we can constructively adjust a x , b l to a = a x r~ x , b = rb 1 , where r G K 
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for O = I and r G XUiT(l + i) for C = J, so as to have M = aOb with a,b eO and 
a primitive (and odd for — J). The argument also shows that once a is adjusted 
to be primitive (and odd if O = J) then b necessarily lies in O. 

Since we are ultimately interested in the similarity submodules, and not so much 
in the actual self-similarities themselves, we have to draw the attention now to the 
symmetries of our maximal orders. The group of units O x of O is, geometrically, a 
finite root system A of typef] -D4 (resp. H4) for J (resp. I). Any self-similarity of O 
that is surjective is necessarily an isometry (norm preserving) and so must map A 
onto itself. Conversely, any isometry which stabilizes A will also stabilize its o-span 
which is O. Thus 

stab o(4) = Aut(A). (28) 

For (9 = 1, Aut(A) is the Weyl group VT(A) of A, which is i? 4 , and consequently 
all elements of Aut(A) + (the orientation preserving part of Aut(A)) are realized by 
mappings M(u,v) with u, v G I x , i.e. units. 

For O = J, [Aut(A) : W(A)] = 6, the additional symmetry being due to the 
diagram automorphisms of D4. In fact, Aut(A) is the Weyl group of F4, and the root 
system of type F 4 can be realized explicitly as 

A := AU{±u±v I u,v G {l,i,j,k}, m/d}, (29) 

i.e. by adjoining to A the elements of J of square length 2. This time, Aut(A) + is 
realized as the set of mappings M(u, v) with either u, v G A or u, v G ( A\ A) / \^2. 
It will be observed that all the elements of A\A lie in the ideal (1 + i) J (in fact they 
are all the generators of this ideal). So, we have proved 

Proposition 2 The orientation preserving self- similarities of O onto itself are pre- 
cisely the maps M(u, v) for u,v G O with \u\ = \v\ = 1 and, in the case O = J, also 
the maps |M(w, v) for u,v G O with |w| 2 = |v| 2 = 2. □ 

We say that an SSM aOb is given in canonical form if a, b G O with a being O- 
primitive (and with a odd if O = J). 

Proposition 3 Similarity submodules a 1 Ob 1 and a 2 Ob 2 written in canonical form 
are equal if and only if both a 1 " 1 a 2 and b 2 b± 1 are units in O. 

6 We are using the symbol D4 both for the root system and for the corresponding root lattice. The 
convex hull of the 24 roots of D4 is the regular 24-cell p. 292], mentioned in an earlier footnote. 
RVM would like to take this opportunity to note that in Q the root system D4 was inexplicably 
left out of the classification of subroot systems of the root system of type H4. All such D4 subroot 
systems are conjugate by the Weyl group of H4. 
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PROOF: Suppose that a 1 C6 1 = a 2 Ob 2 . Then a^ 1 a 2 Ob 2 b^ 1 = O. According to 
Prop. |2], one of two things may happen: 

(i) There is an r G K so that ra^ l a 2 = u E O x (and r~ 1 b 2 b^ 1 = v E O x ). Then, 
ra 2 = a x ii E O gives r E o since a 2 is primitive. Likewise, r~ 1 a 1 = a 2 u 1 E O gives 
r _1 E o, so r G o x C O x and we are done in this case. 

(ii) We are in the case O = J and there is an r G K so that 2ra^ 1 a 2 = (1 + i)u, 
u E O x . This gives 2ra 2 E O and r~ 1 a l = a 2 it _1 (l — i) E O whence 2r, r _1 G Z. 
Thus r = ±1, ±|. Since a x and a 2 are both odd and 2 G (1 + i) 2 O x , none of these 
values of r is possible. 

The reverse direction is clear. □ 
We are now in the position to formulate the main result of this section. 

Theorem 1 The number of similarity submodules of O of a given index is a multi- 
plicative arithmetic function. Its Dirichlet series generating function is given by 

ol ' C K (4s) \l, i/0 = I . ' 

PROOF: As a result of Prop. and Prop. 0, any SSM M. of O can be uniquely written 

as 

M = aOOb, (31) 

i.e. as a product of a right and a left ideal of O, where a is C-primitive (and also odd 
if £> = J). 

Any right ideal dO can be written uniquely as a product c(d)oaO, where c(d) (the 
content of d) is in o and a E O is (9-primitive. In addition, we have the formula 
[O : dO] = [o : c(d)e>] 4 • [0 : aO\. Thus the Dirichlet series for the primitive right 
ideals of O (those aO with a primitive) is the quotient of two zeta functions, 

Co(*)/U4s). (32) 

In the case = 1, the factorization ( pTj ) leads at once to Fj(s) = (Ci( s )) 2 /Cq( t )(4s). 

In the case = J, writing (^) explicitly as an Euler product, we find that the 
contribution for the prime 2 is (l+4~ s ). This corresponds to the fact that the primitive 
right ideals are either of the form a J with a odd or of the form a(l + i)J with a odd. 
The latter ones are to be removed from the counting (because (1 + i)J = JI(1 + i) 
otherwise leads to doubly counting them), and the corresponding Dirichlet series is 
obtained by removing the factor (1 + 4~ s ). This gives the result claimed. □ 
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Results: lattices 



Let us now consider the Hurwitzian case O = J in detail, and also the lattice C = Z 4 . 
If x — (1 + i), we have xO = Ox and also xC = Cx. These lattices are related by 

xO c C C O (33) 

where the integer on top of the inclusion symbol is the corresponding index. 

Lemma 2 If aCb C C, then there exist a l ,b l G (9, a x odd and O * -primitive, 
such that aCb = a 1 Cb 1 . 



Proof: Let aCb C C We can assume, without loss of generality, that a,b G 
and, by using a suitable scaling and the fact that xC = Cx, we may even assume that 
a G O and that a is O-primitive and odd. Since aOxb C aCb C 0, the conditions 
on a already show that xb & O (see Remark 1). Write b = x l c with c G 0. 

Consider [C : a£b] = |a| 4 |b| 4 = ™^ G Z. S ince a is odd, it follows that 4 | |c| 
and hence x | c, because any even element in J is of the form (1 + i) r times an odd 
element P0| . Consequently, b G O and we conclude that aCb C C implies that we 
can rearrange the quaternions in the way claimed. □ 

This provides a link between the SSL problems for £ and for O. The difference 
between the counting arises as follows. The symmetry group of O is isomorphic with 
the Weyl group of F 4 (see above), while that of C, which is the Weyl group of B4, is 
a subgroup of index 3. As we shall show, this only influences the number of SSLs of 
even index when going from O to C. 

Theorem 2 The possible indices of similarity sublattices of hypercubic lattices in J^D 
are precisely the squares of rational integers. The number of SSLs of given index is a 
multiplicative arithmetic function. For the case of J = D\, the corresponding Dirichlet 
series generating function Fj reads 

F , , (CM) 2 (1-2 1 - 2 *) 2 (C(2 g )C(2,-l)) 2 

A ) (l+4" s )C(4s) 1 + 4- C(4s) ' 

The same series also applies to the lattice D4, while for the primitive hypercubic lattice, 
Z 4 , it reads 

FAs) = (1 + J:) • . (35) 

PROOF: The statement about F$(s) follows directly from Theorem [1] and from Eq. (|8|). 
It is rather easy to see from the Euler product representation that precisely all squares 
of integers occur as indices. 
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In order to extend this to £, we have to understand how the different symmetries 
lead to different countings. Assume a£b C £. Due to Lemma |], we may assume that 
a,b G O with a odd and 0-primitive, i.e. we assume canonical form. 

By Prop. [3], aOb = a 1 Ob 1 if and only if there are units u, v G O x with a l = au 
and b Y = vb. However, the unit group of £ is only the quaternion group, Q = 
{±1, ±i, ±j, ±fc}, and Q is a normal subgroup of O x with O x /Q ~ Z/3Z. We may 
take t = (1, 1, 1, l)/2 as a suitable representative of a generator of this cyclic group 
in O x § 26]. Note that, for u,v G O x , we have auCvb = a£b if and only if 
■uv G Q. So the single SSL aOb of may give rise to three different SSLs of £, 
namely to aCb, aCtb, aCt 2 b. Whether or not this happens depends on whether or 
not the latter two are actually in £. 

Now, if a£b and aCtb are both in C, then so is a(C + Ct)b = aOb, whence also 
aCt 2 b C C. Similarly, aCt 2 b C C implies aCtb C £, too. Thus aOb C O gives rise 
to 3 different SSLs of £ if and only if aOb C £. 

However, a(9b C £ implies that [O : a(9b] is even because [O : C] = 2. Conversely, 
[O : aOb] even implies that |a| 4 |6| 4 is divisible by 4, so a or & must be even and 
hence aOb C a;(9 C £. In short, the 3 SSL situation occurs if and only if [O : aOb] 
is even. As a consequence, the counting function for £ is still multiplicative, and the 
modification in the Euler product expansion occurs only in the factor that belongs to 
the prime 2. It is easy to check that the result is that given in the Theorem. □ 

If we take into account that the possible indices are always squares, it is reasonable 
to define the appropriate coefficients as follows, 

m=l 

So, the coefficients actually are 

f c ( m ) = {{£' i s SSL of £ | [£ : £'} = m 2 }\ . (37) 
To simplify explicit formulas here and later on, we introduce the function 

g{n, r) := (r + l)rf + 2 1 ~ (r ± ^ ± T (38) 

for integers r > and n > 1. Note that g>(n, 0) = 1. An explicit expansion of the 
Euler factors now gives the following result. 

Corollary 1 The arithmetic function fj{m) is multiplicative. It is given by 

fM) = I 1 ) \ lf / = 2 AA ■ (39) 
[9\P) r )i if p is an odd prime 
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where r > in all cases. Similarly, f z4 (m) is a multiplicative arithmetic function. It 
is related to /j(m) via 

f ( s f /jM> m od d fAn . 
f z4 (m) = < Q , / \ (40) 
I o • /j(m), m even. 

The first few terms of -Pj(s) read explicitly 

1 n s J L ~ 4 s ~ 9 s ~ 16 s ~ 25 s ~ 36 s ~ 49 s ~ 64 s ~ 81 s ~ 100 s ~ 121 s ~ 144 s ~ V*^ 

while those for F z4 (s) follow easily from (flUl). They are now listed as |29], sequence A 



045771] and [^9], sequence A 035292], respectively. Note that Eq. (|39| ) implies that the 
SSMs of J of index 4 r are unique - they are, in fact, just the 2-sided ideals (1 + i) r JJ. 

Let us now, in line with the previous examples, briefly consider the asymptotic 
behaviour of the coefficients. Since ((s) ^ in {Re(s) > 1}, it is clear that F c (s) is 
meromorphic in the same half-plane, with only one pole which is of second order and 
located at s = 1. This is true both of £ = J and C = Z 4 . Using the Dirichlet series 
(|36l) and applying the results from the Appendix to F c (s/2), we get the following 

Corollary 2 The coefficients f c (m) grow faster than linear on average for large m, 
and we have the asymptotic behaviour 

fc.( m ) - c c' x2 lo gO) (asx ^og) (42) 

m<x 

where the constant is given by 

C C = i»w((.-l)F £ (.)) = ^: J Z1 (43) 

Note that this really is an asymptotic result, and that, numerically, the estimates of 
C c converge rather slowly (from above) to the values given in the Corollary. 

Results: modules 

Let us first state the result for the icosian ring I itself. 

Theorem 3 The possible indices of similarity submodules of the icosian ring I are the 
squares of rational integers that can be represented by the quadratic form x 2 + xy — y 2 . 
The number of SSMs of given index is a multiplicative arithmetic function, its Dirichlet 
series generating function reads 

l[ ) ~ U4s) " C*(4*) (44) 

where K = Q(Y). 
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The proof follows immediately from Theorem [I] in combination with Eq. (|T^). Note 
that the possible indices are just the squares of the possible norms of ideals in Q(r) 
and hence of the form given. Taking this into account, we write 

"M = t^4 («) 

Z — / m zs 

m=l 

in analogy to above, and obtain, by an explicit expansion of the Euler factors, the 
following result (compare [^9], sequence A 035284]). 

Corollary 3 The arithmetic function fi{m) is multiplicative and given by 

g(5,r), ifp = h 

., r . 0, if p = ±2 mod 5 and r is odd 

= ' g(p 2 J), ifp = ±2 mod 5 and r = 2£ ^ ' 

ELo^(P' s )^(P' r_s )' ifp = ±lrnod5 

where always r > and g is the function defined in Eq. $3%). 
The first few terms of F f (s) read explicitly 

F ( e\ - 1 _|_ JiL _|_ JJL _|_ _20_ 4_ 48 , 66 , 80 , 120 , 97 , 120 , 128 , ( A>7\ 

r l\ b ) 1 "T 16 s ~r 25 s 81 s " r 121 s ^ 256 s " t " 361 s 400 s 625 s 841 s ' 961 s " r ' ' ' V '/ 

Let us briefly look at the 10 SSMs of I of index 16. From Eq. ( |T3"D it is obvious that 
they are just the 5 left ideals la and the 5 right ideals al, with suitable generators a 
with N(|a| 2 ) = 16. Note that none of them is 2-sided. 

Finally, we can again determine the asymptotic behaviour along the lines used 
before. F l (s/2) is holomorphic in the half-plane {Re(s) > 2}, with a single second- 
order pole at s = 2. With the results from the Appendix, we then obtain 

Corollary 4 For x — ► oo ; the asymptotic behaviour of the coefficients /j(m) is 

V/i(m) ~ 6(log(r)) 2 ^ 2log ^ ^ 0. 124271 x 2 log(x). (48) 
~^ 5v 5 

At this point, it would be interesting to relate these findings to the corresponding 
ones for the Z[r]-modules C = Z[r] 4 and M. = J[r] — Inf. Since this requires a lot 
more effort than in the previous case (Z 4 versus J), we postpone it, and rather state 
the result for the cubianQ maximal order K in H(Q(y / 2)). 

7 The term "octonian" would be more natural a choice, but it has already been taken! 
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Theorem 4 The possible indices of similarity submodules of the cubian ring IK are the 
squares of rational integers that can be represented by the quadratic form x 2 — 2y 2 . The 
number of SSMs of given index is a multiplicative arithmetic function, its Dirichlet 
series generating function reads 

Fk(s) " "c^fiTT - (49) 

i I 6 i 22 i 32 I 66 I 20 i 192 , 178 , 72 , 120 , 96 , 52 , 
4 s 16 s 49 s ~r 64 s ^ 81 s 196 s " r 256 s 289 s " r 324 s 529 s " r 625 s ^ 

where K = Q(y/2). 

The proof follows directly from the proof of Theorem [I], since literally every step taken 
for the icosian ring translates into one here, with K = Q(a/2), o = Z[a/2], O — K, 
and C = Z[a/2] 4 . Note that now, since 2 is not a prime in Z[a/2] (it actually ramifies 
there), < $ 1 = V2K = (2 + a/2)K is the prime ideal of K sitting on top of the prime 
ideal (2 + \/2)Z[y/2\, and the arguments in the proofs have to be adjusted accordingly. 

One can again work out the coefficients f K (m) explicitly (see |2D|, sequence A 
035285]) 

0(2, r), ifp = 2 

.0, if p = ±3 (8) and r is odd 

j j ^ /y(/) J ./). if p = ±3 (8) and r 2/ 1 j 

Es=o^(P' s )^(P' r_s )' if P = ±1 (8) 



and the asymptotic behaviour, using the Appendix, is 

g(l + 
22^2 



E ^( m ) ~ 15(1 ° g(1 1 ^ ))2 x 2 log(x) ~ 0.374519 x 2 log(x) . (51) 



Concluding remarks 

As we have demonstrated above, the similarity submodules of certain 4D Z-modules 
can be classified by means of algebraic methods based on quaternionic algebras and 
their maximal orders. Together with the results of 0, |6|], this essentially covers the 
cases related to root systems in dimensions d < 4. 

Although we did not emphasize it, one can also determine the actual semigroups 
of self-similarities of these modules explicitly, notably through the canonical represen- 
tation of SSMs (Prop. 0) and their uniqueness up to symmetries (Prop. |3|). We have 
described this in more detail for other cases , and the interested reader will find no 
difficulty to extend that approach to this situation. 
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One application is concerned with the symmetries of coloured versions of the lat- 
tices and modules under consideration. Assume that L has a non-trivial (irreducible) 
point symmetry, and a sublattice L' which is the image of a self-similarity of L of 
index m — [L : L'] > 1. If we assign m different colours to the cosets of L', certain 
subgroups of the point group of L' (which is conjugate to that of L) will give rise to a 
colour symmetry in the sense that their elements induce a unique, global permutation 
of the colours, compare |33[ and references therein. 

This is also closely related to the classification of coincidence site submodules, i.e. 
of submodules that can be written as the intersection of the original module with a 
rotated copy of itself, see for background and some recent results. Here are several 
open questions, particularly in spaces of even dimension, which the above results 
should help to solve for dimension four. 

Finally, one would like to know to what extent a generalization of our results is 
possible. The root lattices seem to form a sufficiently well-behaved class of objects 
to try, and some partial answers on the existence of similarity sublattices and their 
possible indices are given in 15). We are, however, not aware of general results along 
the lines discussed here, i.e. including the determination of the number of SSLs of a 
given index, nor even of a method to overcome the dependence on special features 
such as the arithmetic of quaternions. 
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Appendix 

In what follows, we briefly summarize the results from analytic number theory that we 
need to determine certain asymptotic properties of the coefficients of Dirichlet series 
generating functions. For the general background, we refer to |l| and p6 |. 



Consider a Dirichlet series of the form F(s) = 5^m=i a ( m ) m ~ s - We are mainly 
interested in the quantity A(x) = J2 m <x a ( m ) an d its behaviour for large x. Let us 
give one classical result (based upon Tauberian theorems) for the case that a{m) is 
real and non-negative. 

Theorem 5 Let F(s) be a Dirichlet series with non-negative coefficients which con- 
verges for Re(s) > a > 0. Suppose that F(s) is holomorphic at all points of the line 
{Re(s) = a} except at s = a. Here, when approaching a from the half-plane right of 
it, we assume F(s) to have a singularity of the form F(s) = g(s) + h(s)/(s — a) n+l 
where n is a non-negative integer, and both g(s) and h(s) are holomorphic at s = a. 
Then we have, as x — > oo ; 

A(x) := Va(m) ~ x a (log(x)) n . (52) 

a ■ nl 

m<x 

The proof follows easily from Delange's theorem, e.g. by taking q = and u = n in 
Tenenbaum's formulation of it, see 0, ch. II. 7, Thm. 15] and references given there. 



Note that Delange's theorem is a lot more general in that is still gives results for 
other local behaviour of F(s) in the neighbourhood of s — a, in particular for n not 
an integer and even for combinations with logarithmic singularities. Let us also point 
out that there are various extensions to Dirichlet series with complex coefficients, e.g. 



Thm. 1 on p. 311 of [24]], and even stronger results (with good error estimates) for 
multiplicative arithmetic functions a(m) with values in the unit disc, see 0, ch. I, § 
3.8 and ch. Ill, § 4.3] for details. 

The critical assumption in Theorem |5] is the behaviour of F(s) along the entire line 
{Re(s) = a}. In all cases that appear in this article, this can be checked explicitly. 
To do so, we have to know a few properties of the Riemann zeta function, ((s), and 
of the Dedekind zeta functions of Q(r) and Q(y/2). It is well known that ((s) is a 
meromorphic function in the complex plane, and that is has a sole simple pole at 
s — 1 with residue 1 |l], Thm. 12.5(a)]. It has no zeros in the half-plane {Re(s) > 1} 
|34] , ch. II. 3, Thm. 9]. The values of £(s) at positive even integers are known Jl|, Thm. 
12.17] and we have 



n 2 tt" 



C(2) = ¥ , C(4) = - (53) 
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This is all we need to know for this case. 

The Dedekind zeta function of K = Q(r) has some similarly nice properties. It 
follows from |I5|, Thm. 4.3] or from |3(| §11, Eq. (10)] that it can be written as 



Cq (t )(*) = C(s)-L(s, X ) (54) 

where L(s, x) is the L-series of the primitive Dirichlet character [l|, Ch. 6.8] \ defined 
by 

f 0, n = (5) 
X (n) = I 1, n = ±1 (5) (55) 
(-1, n = ±2(5) . 

Since x is not the principal character, L(s,x) = Ylm=i x( m ) m ~ s is an entire function 
|l], Thm. 12.5]. Consequently, Cq( t )( s ) is meromorphic, and its only pole is simple and 
located at s — 1. The residue is L(l, %) and from |35|, Thm. 4.9] we get 



2 Mr) 
v/5 



res s=1 ( Q(T) (s) = L{1, X ) = %^ ^ 0.430409. (56) 



Since L(s,x) ^ for Re(s) > 1, see |3^, p. 31], Cq( t )( s ) cannot vanish there either. 
Also, one can again calculate the values of Cq(t)( s ) a ^ positive even integers. This is 
done by means of the functional equation of L(s, x) P- 30] and the knowledge 



of the values of L-functions at negative integers in terms of generalized Bernoulli 
numbers [^5], Thm. 4.2], see |35|, Prop. 4.1] for a formula for them. Working this out 



explictly for s = 2 and s = 4 gives 

27T 4 4-7T 8 

Cqm(2) = 757E ' Cqm(4) = TeWf- (57) 

Let us add that these results, and those to follow, can also be found, in rather explicit 
form, in §9 and §11 of |3"6f . 

In the same way, one can determine the zeta function of Q(\/2) and its properties. 
One has Cq^)^) = C( s )L(s, x)i now with the primitive Dirichlet character x{ n ) = 
0, 1, —1 for n even, n = ±1 (8), n = ±3 (8), respectively. This zeta function has again 
only one simple pole, at s = 1, with residue L(l,x) = log(l + V2)/\/2 ~ 0.623225. 
Finally, we have 



W 2 ) - . C Q( v5)(4) - gg^^- (58) 
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